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Abstract 


Identifying insecticide resistance mechanisms is paramount for pest insect control, 
as the understandings that underpin insect control strategies must provide ways 
of detecting and managing resistance, thus study identified the ability of 
Carboxylesterase (CarEs) genes in conferring resistance to commonly used 
Article Info pesticide by examining the expression of the gene in the susceptible and resistant 
strain of Zonocerus variegatus, after Cypermethrin bioassay. Adult Z. variegatus 


Volume 3, Issue 4, October 2021 | as sampled from Savannah area of Oyo State, Nigeria, samples were collected 


Received : 03 April 2021 and preserved in “RNA Later” for RNA extraction. Thetwo groups of Z. variegatus 
Accepted : 14 September 2021 were later subjected to quantitative quantitative Polymerase Chain Reaction 
Published: 05 October 2021 (qPCR) using two primers. CarEs gene expressions have both negative and 


doi: 10.33472/AF]BS.3.4.2021.76-87 positive AAct value when LmcesA 20 primer was used. However, no significant 
difference in expression level was observed between groups for both primers 
LmcesE9(F =0.84, p >0.05) and LmcesA 20(F =0.499, p >0.05). Thestudy concluded 
that since the expression level observed for CarEs gene was the same using 
LmcesE9 in both groups, LmcesE9 may not function for detoxifying Cypermethrin 
In Z. variegatus. LmcesA 20 showed a little variation in expression level in both 
strains having higher fold change in resistant strain than the susceptible strain, 
thus these members of CarEs genes need to be further verified for their specific 
role in conferring resistance to use of pesticides. 
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1. Introduction 


Carboxylesterases (CarEs) play significant roles in metabolism of certain hormones and detoxification of 
dietary and environmental xenobiotics in insects. These classes of enzymes are responsiblefor hydrolysis of 
chemicals which contain a functional group of carboxylic acid ester, amide, and thioester (VWheelock ef al., 
2005; Yang, 2012). Many insect CarEs areassociated with insecticideresistance or hormone and seniochenical 
metabolism (Oakeshott et al., 2010; Oakeshott et al., 2013). They are involved in the detoxification of many 
insecticides, including synthetic pyrethroids (SPs), carbamates (CBs) and organophosphates (OPs) which 
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havebeen shown to beinvolved in development of resistanceto insecticidegroups SPs, CBs and OPs (Cui & al., 
2011; and Wu e al., 2011). 


Zonocerus variegatus, (Orthoptera: Pyrgomorphidae) genus of grasshoppers (Caelifera) in the family 
Pyrgomorphidae (Roskov é& al., 2011). It has been regarded as the main pest of food crops such as citrus, 
banana, Cassava, cocoa, coffee, cowpea and yams in West and Central Africa’s moist lowland forest and 
savannas (Modder, 1994; and Kekeunou & al., 2006). Over the last 20 years, Z. variegatus outbreak have 
increased dueto the high production rate of cassava and due to widespread or distribution of Chromolaena 
odorata (L.) (Bamideleand M use, 2012). By thetimethefirst, second and third instars liveon weeds such asC. 
odorata and A spilia africana (Asteraceae family), thefifth and sixth instarsfeeds on M . esculenta (Euphorbiaceae) 
(Chapman et al., 1986). 


Cypermethrin belongs to a class of insecticide known as synthetic pyrethroids. It has a high level of 
insecticidal activity, low effect on class Aves and Mammalia. Itis most preferableto other form of insecticides 
because it is less toxicto non target organisms and since others pesticides have been gradually phased out. It 
belongs to aclass of insecticidethatis widely used by farmersto control Z. variegatus, and it can also haveeffect 
on non-target organism like human being. 


CarEs pyrethroid metabolism has been the subject of several studies due to its widespread usage in 
agriculture (Hodgson and Levi, 1996) and public health (Takken, 2002). CarEs plays an important role in 
detoxifying pyrethroids, which was demonstrated practically in recent rodent studies (Gaughan et al., 1980). 


CarEs play akey rolein the metabolism of xenobiotics including insecticides in insects. The understanding 
of expression patterns of such detoxifying gene and as well as it the effect of insecticides on its enzyme 
activities are necessary to clarify thefunction of this gneredevant to insecticides-detoxifying process. But little 
information is availablein the variegated grasshopper Z onocerus variegatus (L.) 


Insects have developed resistance to almost all classes of chemical pesticides through detoxification 
mechanisms involving anumber of detoxifying enzymes, including glutathione S-transferases (GSTs) (Ranson 
et al., 2001). B-esterases, including Cholinesterases and Carboxyesterases, are typically inhibited by OPsdue 
to the extremely slow dephosphorylation of tetrahedral intermediates produced between OPs and aserine 
residue at their activesites (Fukuto, 1990). Thesethreetypes of enzyme esterases (through ester hydrolysis), 
cytochrome P450 monoxygenases (through oxidation), and glutathionetransfarases (through ester hydrolysis) 
are primarily used to turn insecticide into less harmful products (H ol!ingworth and Dong, 2008). 


Pyrethroids and Carbamates are believed to be hydrolyzed by CarEs. The quantitative approach Is 
determined by the overproduction of CarE proteins by gene amplification and transcriptional up-regulation. 


2. Materials and methods 


2.1.1 nsect collection 


Z.variegatus weresampled in Savannah area of Oyo state using a sweep net (Kemp ef al., 1990; and Cigliano 
et al., 2000). They weretransferred to insect cage of dimension 29.5cm x27.5cm x40.5cm for grouping. They 
werefed with cassava leaves for the period of collection before exposureto Cypermethrin. 


Thespecimen collected weretransported to the Entomology laboratory, wherethey weresorted In groups 
into hoarding plastic cage of dimension 10cm x12.5cm x6cm for Cypermethrin bioassay. 


2.2. Exposure to chemical and RNA extraction 


Groups of adult Z. variegatus weretested for insecticideresistivity and some group without insecticide are 
used as the control. Z. variegatus and the cassava leaves were sprayed with different concentration of the 
recommended treatments (Cypermethrin) in different groups except for the control and observed for 1h. The 
knockdown timeof individual grasshopper were recorded in minutes, thedead ones arecollected after 1h of 
exposure and deemed insecticidesusceptibleif after mechanical stimulation, they remain motionless. The Z. 
variegatus alive after one (1) hour were scored as cypermethrin-resistant. The resistant and susceptible Z. 
variegatus arepreserved in “RNAlater” (Qiagen) for RNA extraction. 


Total RNA was extracted from three samples of susceptible and resistant strains of each concentration of 
chemicals used using Trizol reagent (Takara, Dalian, China) and treated with RNase-free DN ase! (Promega 
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Corporation, Madison, WI), purity of the RNA was quantified by measuring UV absorption using a 
spectrophotometer, absorbance is measured at 260 nm (A260) and 280 nm. The first-strand CDNA was 
synthesized by reverse transcribing from 1.5 wg of total RNA using MLV reversetranscriptase (Takara). The 
relative expression levels of 25 CarE-likegenes of thelocust by real-time quantitative PCR (qPCR) with Z-actin 
act as areferencegene. 


2.3. Primer selection 


Two primers were used for this study LmcesA 20 and LmcesE9 based on specific band showed by gel 
electrophoresis.PCR amplicon electrophoresis was carried out by size fraction on 2.0% Agarose gels. 
Electrophoresis was used to determinethe quality and integrity of the RNA and thebinding of primers with 
CDNA by fractionation on 2.0% agarose gels. Agarose gels were prepared by dissolving and boiling 0.7 g 
agarose in 35 ml 0.5 x TBE buffer solutions. Electrophoresis was done for 75 °C for 30 min. The primer 
sequences and the expected sizeof each PCR product areshown in Table 1. 


Table 1: Table of primer selected to test for Carboxylesterase gene in Z. variegatus 


Gene name Application of primers Sequences (59-39) Product size (base point) 


[earenecrerimmerece 


B-actin F: CGAAGCACAGTCAAAGAGAGGTA a 


Lmces E9 F: CAGAACCTCCTGTTGGAACACA 


R: CAGAGCATCTCTTACACCATTCCAT 





2.4.Q uantitative- PCR verification of CarEs genes associated with C ypermethrin resistance 


TheqPCR was performed by using threebiological samplesin each groups, Relative expression was calculated 
using the 2“4*t method. Analysis of Variance was performed to determine significance of difference in the 
relative expression of thetwo (2) CarEs genes between the susceptible and resistant strain. 


3. Results 


3.1. Primers selection 


Electrophoresis showed B-actin to besignificantly expressed than the other primers and two primersLmcesA 20 
and LmcesE9 were chosen from five primers selected fromL. migratoria CarEs genes based on specific band 
showed. The expression levels of the primers se ected werevery low inthesedected sample. 
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Table 2: RNA concentration measured by spectrophotometer for each of the samples extracted 


} Sample 1D glut | A280 } 260/280 260/230 ) Const. [Cursor Pos. Pos. | Cursor [Cursor abs.| 340 raw 


ome 48 175. | 175.32 


3 [ia | 0 
rae [am [ee [oo fom] =| = | foo 
roves [wn som [awe [ow [ae] @ | = | oo [ow 


|e feefaefoe foes = poe Lon 


3.2. Cycles threshold (Ct) 


Cyclethreshold (Ct) values for each of the samples were calculated according to the values generated from the 
analysis. All the Ct values were subjected to One Way Analysis of Variance (ANOVA), results showed 
that there was no significant difference between the susceptible strains and resistant strains for LmcesE9 
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(F =0.734, p >0.05) and LmcesA 20 (F =0.730, p >0.05). There was no significance difference between both 
strains using LmcesA 20 (F =0.509, p >0.05) and LmcesE9 (F =0.281, p >0.05). Amplication curve for all the 
samples using the ct valuefor both primers were plotted and shown in Figures 1 and 2 respectively. 


_ Plateau Phase 


Ysampie 


SS SS SS SSS SS SOS TS STS BS SSeS co 


Initiation phase .@ 


0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 4042 
Cycle 


Figure 1: Amplification curves for LmcesA 20. This graph shows the fluorescence increase for different 
samples with the number of cycles. The threshold was set above the limit of detection, but well below the 


stage where the amplification rate slows down 


C. value 


Sample 


Fluorescence 


Initiation phase” 
ry, 


<u NTC 


0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 4042 
Cycle 


Figure 2: Amplification curves for LmcesE9. This graph shows the increase of fluorescence with the number 
of cycles for different samples. The threshold was set above the detection limit but well below the plateau 
phase where the amplification rate slows down 


Note: NTC means Null Template Control 
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3.3. Comparative quantification using AACt (24) 


In using LmcesE9the highest AAct value was obtained from susceptible strain DO.9 (-22.30) and lowest value 
from A 0.48 (-14.18) from resistant strain, likewise using LmcesE9 highest AAct value was from A 2.9 (-6.21) 
from resistant strain and lowest AAct valuewas from D 0.9 (-22.3) from resistant strain as well. Ct valuevaried 
significantly in both the primers, meanwhile AAct value of 0.25 indicated 4X decrease in transcript level 
compared to housekeeping gne. Oneway ANOVA wasused to compared significancedifference between the 
AAct value in susceptible strain and resistant strain using both primers, values gotten showed no significant 
differencein susceptiblestrain and resistant strain for both primers LmcesE9 (F =0.84, p >0.05) and LmcesA 20 
(F =0.499, p >0.05). Fold change expression using LmcesE9 range from 5163793-33689 for susceptiblestrains 
and range was 1100707-18561 for resistant strains. For LmcesA 20 fold change range from 0.54-1.71 for 
susceptible strains and 1.69E-7 - 74.03 for susceptible strains. 


The relative expression of CarEs genes measured by qPCR was represented by diagram of fold gene 
expression. It was plotted using the descriptive statistics calculated, the valueof calibrator used was obtained 
by taking the average of all control values, for LmcesE9 it was 12.68. The average value of 2““‘for all the 
samples in both treated (Susceptibleand resistant strains) and control group were calculated for LmcesE9, it 
was 498.60 and 0.45 respectively, whilethe standard deviation in the treated was 345.17 and control was 0.77. 
Thestandard error calculated was 199.28 in thetreated group and 0.45in thecontrol group. Theplot of thefold 
gene expression for LmcesE9 shows high fold gene expression level in the treated compared to the control 
group (Figure3). Thedescriptivestatistics calculated for LmcesA 20 showed the valueof calibrator used and it 


Table 3: Calculated 2:“4‘' for LmcesE9 primers with respect to the gene of interest and housekeeping gene (B-actin1) 


E9 CT1| E9 CT2 HKG HKG HKG AACt 
eat ct Avg. Ct a 


ro | | om [ow [ww [ww [ow [ 
ona | aw | ew [aw | vm | se | arm | am | ive | woesien 
rose | am | aman [nw | mm | mw | am | oe en 


A/ 0.48 11.49 11.49 11.49 12.99 12.99 12.99 ces -14.18 18561.17 


A/ 0.9 


fe taste test eal et ae fo ee 
moa [ en [on [en [uw [ese [om [oy [oes | mmm 
vo [re [7 [oe | om [om [am | sam] ene 
ae [es [sf em [we [sw [ao [oe [so 


Note: GOI means gene of interest, HKG means house keeping gene. 
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Table 4: Descriptive analysis for AAct values of LM CESE9 primers in each sample and average value of 


calibrator generated from control group (normalizer) 


Calibrator (average ACt control group) a 
Descriptive statistics Sy 
vem | Ame 


a 


Note: SD mean Standard Deviation; and SE means Standard Error. 


Table 5: Calculated 2-“4CT for LmcesA 20 primers with respect to the gene of interest and housekeeping gene (B-actin2) 


A20 A20 A20 HKG HKG HKG ACt AACt 2*- 
a aa Ct Ct2 Avg. Ct (AACt) 


renas | wm [om fam [we | mw | we | ow | om | eo 
Temas [pa [aw [ow [vs [vs [we [ao | om [ae 
Foreae [me [nm [mm [wo [wo [uo | am | oe [om 
fos | en fon|en [oe | oe] ow | am | on | on 
pene [om [om [om fw [or [er [ee [a [oe 
fess [on [nn [ae [ew [ew [ew [oe [om [om 
wes [sn fon [sm [aw [ow | om [op [sm [on 
rons | ew [ow [om [as | asus [on | on | om 
pene [om [om [oe [an [om [em [ow [om om 
fens [sue [ase [nae [oe [ao [ae [nw [2s [soe 
pans [ea [ew [aa [se [se Pose [aor [a [ ow 

pass [ew [om [om es [ae [nm [oes [on [em 
pre | mo [aor] aoe [awe | aw | ow | ae | oe | ow 


was obtained by taking the average of all control values as well, for LmcesA 20 it was -4.04, the averagevalue 
of 2““'for all thesamples in both treated (susceptible and resistant strains) and control group was observed 
to be 19.72 and 918.87 respectively, while the standard deviation in the treated was 9.28 and control was 
1577.23. The standard error calculated was 5.36in thetreated group and 910.61 in thecontrol group. Theplot 
of thefold geneexpression for LmcesA 20 shows high fold gene expression level in the control group compared 
to thetreated group (Figure 4). 
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Table 6: Descriptive analysis for AAct value of LmcesA 20 in each sample average value of calibrator generated 


from control group (normalizer) 


caer acengen [a 
ae as 
eee 


Note: SD mean Standard Deviation; and SE means Standard Error. 


®/ control, 55.36 


S 
2 
“ 
“4 
w 
— 
2. 
x 
@ 
w 
— ’ 
@ 
bo 
a) 
5 
as 


® treated, 0.33 


Groups 1.00 


Figure 3: Fold gene expression level for control and treated group using LmcesA 20 
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Figure 4: Fold gene expression level for control and treated groups using LmcesE9 
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3.4. Expression profiles of CarEs genes 


Theexpression profiles of CarEs genes measured by gPCR in all the evaluated samples were represented by a 
hierarchical clustering in both primers. It showed therdationship between expression of CarEs genes In both 
susceptible and resistant strains for all the evaluated samples. Using LmcesE9 primers the expression leve for 
CarEs genes were viewed to beat thesame level for both the susceptiblestains and resistant strains, as well as 
the control group but theonly out group which the expression level was different was said to be NTC (No 15) 
inwhichthetemplatecDNA sampleis absent, but only provides a measure of reaction set up contamination 
(Figure6). Using LmcesA 20 primers expression level of CarEs genes werealso observed to beat thesamelevea 
for both susceptible and resistant strains but differences were observed in the control group (No 1) and 
susceptiblestrainA 2.9(N 014), Figure5. 












































Figure 5: Dendrogram representing expression profiles that shows the relationship of CarE genes observed 
In both susceptible and resistant strains using LmcesA 20 primers measured by qPCR 


Figure 6: Dendrogram representing expression profiles that shows the relationship of CarE genes observed 


in both susceptible and resistant strains in Z. variegatus using LmcesE9 primers measured by qPCR 
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4. Discussion 


4.1. | dentification of CarEs genes 


ThecDNA of Z. variegatus subjected to gPCR using primers from L. migratoria indicated expression of two 
CarEs alone, though CarEs gene maybe more than two as it was revealed from literature review of other 
organism likeL. migartoria in which about twenty five CarEs genes werebeing identified (Zhang et al., 2013). 


N egative value for AAct shows there is an up-regulation of genes, all samples with negative AAct values 
have an up-regulation of CarEs genes in them indicating the CarEs genes is well expressed/ amplified in 
LmcesE9 primers, but results showed using LmcesA 20 that some samples have the CarEs genes being up- 
regulated whilesome being down-regulated. Negative sign for AAct givevalue greater than 1 for fold change 
expression which automatically indicateup-regulation of thegene, whilepositivesign for AAct giveafraction, 
|.e.,avaluelessthan 1 for fold changeexpression which indicate down-regulation. 


Expression leved of LmceE9 was revealed to behigh in both susceptibleand resistant strains of L. migratoria 
intheair responseto malathion, (Zhang & al., 2010). This study also revealed thesame expression pattern using 
LmcesE9in their resoonseto Cypermethrin. The expression level observed in LmcesE9 was said to behighin 
the susceptible strains than theresistant strains. For expression level using LmcesA 20 it was observed higher 
expression was obtained in the resistant strains than the susceptible strains. LmcesE9 was found to behighly 
expressed in both the susceptible and resistant strain and susceptible strain was significantly up-regulated 
with about 5163793 fold, but in LmcesA 20 CarEs identified was relatively low in expression though the 
strains differ in their level of expression, highest fold was obtained in the resistant strain with 74.3 fold. 
Among the two CarEs genes LmcesA 20 tends to show variation between thetwo strains and will be good for 
further studies on Z. variegatus resistanceto pesticides. 


Increased production of CarEs due to gene amplification or gene transcription up-regulation has been 
found to confer insecticideresistancethrough increased hydrolysis and/ or insecticidesequestration in resistant 
insects (Smal! and H emingway, 2000; and Cui & al., 2007). This study observed expression level of CarEs gene 
but not over expressed as compare to other studies. M ost of the previous work compared changes in enzyme 
activity between susceptible and resistant strains. Reyes e al. (2011) demonstrated that higher activity inthe 
responsive strain was observed than in the codling moth’s OP-resistant strains which support that same 
degree of expression may not besignificantly different in susceptible and resistant strains. 


In this study, both theCarEs gene identified had almost thesameleve of expression in both thesusceptible 
and resistant strains though the gene was highly expressed using LmcesE9 primers, which signifies variation 
in the CarEs genes of insects and pointing to other functions than resistivity to Cypermethrin. Therelatively 
expression of CarEs gene using LmcesE9 and sameleved of expression in both strains may function in resistance 
to other pesticides other than Cypermethrin (cross resistance). 


4.2. Expression profiles of CarEs genes in Z. variegatus 


This study verified the activities of CarEs genes that may be associated with Cypermethrin resistance and this 
was aS aresult of recent reports on the resistivity of Some Z. variegatus to insecticide, and the fact that these 
genes have been earlier reported in other insects that they confer resistance. Wherefore it Is often not clear 
whether these genes Is expressed in Z. variegatus and the role it played in conferring resistance to use of 
pesticides, thus making theselectivity and sensitivity of our GPCR method, and Cypermethrin as pesticideas 
itis bang commonly used, ‘especially owing to the fact that these organisms are causing destruction to farm 
produces which its effects need to be curtailed. 


Expression leva for CarEs gene was observed in Z. variegatus for both susceptible and resistant strains. The 
expression of CarEs genes in Z. variegatus suggested some potential roles itis playing in insects. Therefore, the 
expression pattern in both susceptible and resistant strains and Its distribution givea basic knowledge that all 
the CarEs have different expression pattern and theroleeach Is playing differ from each other in Z. variegatus. 
Expression profiles areimportant for the purposeof transcription profiling and Class discovery of CarEs. The 
two CarEs gene expressed was probably not of the same class because different expression leved. 
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4.3. Comparing expression profiles in both susceptible strain and resistant strain 


This present study focus on the quantitative expression profiles of CarEs genes in susceptible strain and 
resistant strain of Z. variegatus; the first finding, examined the expression of CarEs using different primers 
which indeed varied significantly while the second finding examined the expression profiles for both the 
susceptible strain and the resistant strain in all the Cypermethrin concentration used, which indicate no 
significantly difference. Relative higher expression level was obtained using LmcesE9 primers but lowly 
expressed using LmcesA 20 this indicate specificity in the use of primers. The tissue in which the analysis 
compared the expression pattern may seem to contribute to the expression pattern in thestrains because of 
assumed diversefunctions of CarEs genes. 


5. Conclusion 


Thus, the study concluded only LmcesA 20 have potential rolein conferring resistanceto use of Cypermethrin 
in Z. variegatus becauseof different expression level in susceptible and resistant strains which can possibly be 
agood candidate of the detoxification genes. LmcesE9 havesame level of expression in both theresistant and 
susceptible strains and so it point to other biological roles performed by CarEs genes, e.g., production of 
hormones and degradation of pheromones. 


6. Recommendations 


From this study, the following aretherefore recommended in order to obtain data that can be used to verify the 
clades of genes involved in conferring resistance: 


1. Sequencing of CarE genes in Z. variegatus, Study should comparedifferent expression level of CarEs genes 
in other tissues other than the Limb. 


2. Clarify the CarE activity variations correlated with genetranscript variations, moreCarE genes should be 
further investigated. 
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